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Abstract. We evaluate the effects that variations in He content have on bolometric corrections and r eff -colour relations. To 
this aim, we compute ATLAS9 model atmospheres and spectral energy distributions for effective temperatures ranging from 
3500 K to 40000 K for dwarfs and from 3500 K to 8000 K for giants, considering both "He-non enhanced" and "He-enhanced" 
compositions. The considered variations in He content are of AY = +0.1 and +0.2 for the metallicity [M/H] = +0.5 and 
AK = +0.1 for [M/H] = -0.5 and —1.5. Then, synthetic photometry is performed in the UBVRIJH K system. We conclude that 
the changes in bolometric corrections, caused by the adopted He-enhancements are in general too small (less than 0.01 mag), 
for both dwarfs and giants, to be affecting present-day tables of bolometric corrections at a significant level. The only possible 
exceptions are found for the l/-band at T s g between 4000 K and 8000 K, where \ABCu\ amounts to ~ 0.02 mag, and for T e g 
equal to 3500 K, where |ABCs,| values become clearly much higher (up to 0.06 mag for passbands from U to V). However, 
even in the latter case the overall uncertainty caused by variations in the He content may be not so significant, because the 
ATLAS9 results are still approximative at their lowest temperature limit. 



1. Introduction 

Over the last years, many different problems have prompted the 
computation of stellar evolutionary tracks for different values 
of initial He content. Just to mention a few of them: 

1) Evidence has been found for significant variations in He 
content in some globular cluster like a> Cen, NGC 2808 and 
M 13 (Piotto et al. 2005; Lee et al. 2005; D'Antona et al. 2005; 
Caloi & D'Antona 2005). These variations suggest that the re- 
lationship between He and metal content was not univocal dur- 
ing the first period of chemical enrichemnt of the universe. A 
univocal relationship, instead, has so far been assumed in most 
grids of stellar models applied to the study of old stellar popu- 
lations. 

2) The primordial He content has been recently revised up- 
ward, from Y p ~ 0.235 to Y p = 0.248 + 0.001, after the WMAP 
mission (Spergel et al. 2003, 2006). Many grids of models for 
population II stars have been computed for Y values lower than 
the WMAP one. 

3) On one side the helium content in five Hyades binary 
systems (Y = 0.255) is lower than expected from their super- 
solar metallicity, pointing to a value dY/dZ of the order of 1 
(Lebreton et al, 2001), whereas other observations either indi- 
cate higher values, dY/dZ ~ 2 — 2.5 (e.g Jimenez et al. 2003 
from K dwarf stars in the Hipparcos catalog; Peimbert et al. 
2002 from extragalactic HII regions), or fail to constrain it to 
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a significant level (dY/dZ =3 + 2, Pagel & Portinari 1998). 
Grids of stellar models for population synthesis (e.g. Bertelli et 
al. 1994; Girardi et al 2000), instead, in general use high dY/dZ 
values in order to fit both the primordial and the solar initial He 
content. 

These aspects have prompted us to start a large project for 
the computation of stellar tracks covering a large region of the 

Y - Z plane (Bertelli et al. in preparation). Once ready, these 
tracks will allow us to model stellar populations at any inter- 
mediate Y, thus taking into account the changes in lifetimes, 
luminosities and r e jf that follow from a varying Y. 

However, before stellar evolutionary tracks and isochrones 
are compared with observations, they have to be converted to 
magnitudes and colours via bolometric corrections (BC) and 
colour- r e ft relations. The latter may also be affected by the 
changes in He content, and the purpose of this paper is exactly 
to evaluate how much. To do so, we first compute energy dis- 
tributions for a few selected chemical mixtures with different 

Y (Sect. |2), and then perform synthetic photometry on them 
(Sect. [3). The results, in terms of changes in BCs and colours, 
are discussed in Sect. [4] 

2. Synthetic spectra for He-enhanced 
compositions 

Small grids of ATLAS9 model atmospheres and energy distri- 
butions (Castelli & Kurucz 2003) were generated for different 
sets of metallicities and enhanced helium contents. For con- 
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sistency reasons between continuous and line opacities, new 
opacity distribution functions (ODFs) were computed for each 
chemical composition having enhanced helium abundance. The 
DFSYNTHE code (Kurucz 2005; Castelli 2005) was used to 
this purpose. 

The solar and scaled-solar abundances selected for this 
study are based on the solar chemical composition from 
Grevesse & Sauval (1998). They are the same ones used by 
Castelli & Kurucz (2003) for the ODFNEW grids of models 
and fluxe^]. In terms of fractional mass, the abundances are 
X = 0.735, Y = 0.248, Z = 0.0170 for the solar case. The solar 
and scaled-solar abundances will hereafter be mentioned as the 
AT = case. Then, for 3 different values of metal content, we 
have computed energy distributions for the following mixtures: 

- for [M/H] = -1.5: AT = and AT = 0.1. 

- for [M/H] = -0.5: AT = and AT = 0.1. 

- for [M/H] = +0.5: AT = 0, AT = 0.1 and AT = 0.2. 

The [M/H] = -1.5 and [M/H] = -0.5 spectra aim to probe 
the effect of He at globular cluster metallicities, whereas the 
[M/H] = +0.5 ones serve to probe the potential effect at the su- 
persolar metallicities found in giant ellipticals, for which mea- 
surements of the He content do not exist. The effects at solar 
metallicities are of course derivable by interpolation between 
the [M/H] = -0.5 and [M/H] = +0.5 cases. 

Then, for each one of these chemical mixtures, we com- 
pute energy distributions for a sequence of dwarfs and giants at 
several T e ff values. They are: 

- Dwarfs: with logg = 4.5, and for T eff =3500, 4000, 5000, 
6000, 8000, 12000, 20000, and 40000 K. 

- Giants: with logg = 1.5, and for r eff =3500, 4000, 5000, 
6000, and 8000 K. 

As an example, Figure Q] compares spectral energy distri- 
butions, differing only for the He content, for a relatively cool 
dwarf of intermediate metallicity. In the top panel, the upper 
continous lines indicate the emergent flux due to the only con- 
tinous opacities, while the lower lines are the emergent flux due 
to both continous and line opacities. The He-enhancement has 
a modest impact on the emergent spectra. This is evident in the 
bottom panel of Fig. Q] where the quantity 1 - F & Y = 01 /f^ y=0 
is plotted. The differences between the two spectra amount to 
just a few percent, which translate in maximum changes of just 
a few hundredths of magnitude in bolometric corrections (see 
Sect.[3]below). 

Moreover, some of the differences seen in the bottom panel 
of Fig. [T]are of no concern because they appear at spectral re- 
gions where the emergent flux is very small (for instance, for 
A < 400 nm in the figure). In order to better illustrate the differ- 
ences in the computed spectra which are due only to the varia- 
tion in the He content, Fig.|2]presents a complete series of plots 
of the quantity 6F A , defined as 



5F> = 



r A 



pAy=o.i 

r A 



7 AY=0 



(1) 



1 The ODFNEW spectral energy distributions from Castelli & 
Kurucz (2003), as well as the He-enhanced ones presented in this pa- 
per, are available at 

http : //wwwuser . oat . ts . astro . it/castelli/grids . html 




Fig. 1. Top panel: Spectral energy distribution for a star with 
T eff = 3500 K, logg = 4.5, and [M/H] = -0.5, for both the 
AT = and AT = 0.1 cases (full and dotted lines, respectively). 
The two upper lines compare the continua, the two lower 
lines compare the emergent fluxes (f& Y=Ql and F^ Y=0 ). Bottom 
panel: the relative difference between the above spectral en- 
ergy distributions, illustrated by means of 1 - F & Y=0 - 1 / f& Y =° _ 
The plot also indicates the approximate location of Johnson- 
Cousins UBVRIJHK pass-bands. 



where F^° is the maximum flux of the F^ Y=0 spectrum. By 
plotting the quantity SFa, we evidence only the differences that 
occur in the spectral region which is more relevant in terms 
of flux. This allows a quick evaluation of the changes that are 
potentially more important to the photometry. Of course, dif- 
ferences between the AT = and AT > cases occur over the 
complete range of A. 

3. Synthetic photometry and results 

We have performed synthetic photometry for the above- 
mentioned energy distributions using the same formalism as 
in Bessell et al. (1998) and Girardi et al. (2002). Since we are 
just interested in the changes that the enhanced He can have in 
the synthetic photometry, the equation to be used is: 



ABC, 



Sa 



-2.5 log 



f=°S A dA 



(2) 



where Sa is the total throughput in the filter under considera- 
tion, defined in the interval [Ai , A2] - These ABCs t directly tell 
us the effect of He-enhancement on the absolute magnitudes. 
The effect on colours can be simply derived by the differences 
in ABCs 1 f° r two filters. 

Figures[3]to|6]illustrate the behaviour of ABCs A as a func- 
tion of T e ff, for both dwarfs and giants, for all [M/H] and AT 
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Fig. 3. The ABC Sa quantities as a function of r eff , for both Fig. 5. The same as Fig. [3] but for [M/H] = +0.5 and AY = 0.1. 
dwarfs and giants of metallicity [M/H] = -1.5 and AY = 0.1, 



and for some of the Johnson-Cousins-Glass UBVRIJHK fil- 
ters. The small crosses are the ABCs , values effectively com- 
puted in this work; they are linked by natural spline curves just 
for the sake of a better distinction between the different filters. 
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Fig. 4. The same as Fig. but for [M/H] = -0.5 and AT = 0.1. 



values considered in this work, and for the specific case of 
Johnson-Cousins-Glass UBVRIJHK filters. The filter curves 
were taken from Bessell (1990) and Bessell & Brett (1988). 
The same data are tabulated in Table 1, and is provided in elec- 
tronic form at http://pleiadi.oapd.inaf. it 



T„ (10= K) 

Fig. 6. The same as Fig.EJ but for [M/H] = +0.5 and AY = 0.2. 

4. Discussion and conclusions 

As can be readily seen in Figs. [3] to [6] the effect of He- 
enhancement in the BC is overall quite modest. The most re- 
markable result for all cases considered here is that | ABCs J 
are smaller than 0.009 mag for all stars with T e ff > 5000 K and 
for all pass-bands redder than U. The typical \ABCs , I values 
in these cases are even smaller, of the order of 0.005 mag. In 
general the corresponding shifts in absolute magnitude are well 
below the typical errors in photometric observations. Since in 
most cases the ABC$ , behave in a similar way for different fil- 
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ters, the effects in colours are of even smaller magnitude. It is 
clear to us that the effect of He-enhancement is small enough 
to be neglected in these cases. 

Significant values of IABC5J are met just in a few situa- 
tions; namely for the U filter and at intermediate values of 
i.e. between 4000 and 8000 K, \ABC Sl \ become slightly higher 
but anyway still of the order of 0.02 mag. This is already an 
effect that could be detected in the (as far as we know, rare) 
case of high-precision photometry in the U passband. The con- 
sidered range of T e ff is high enough to include the turn-off re- 
gion of metal-poor globular clusters, and part of their horizon- 
tal branch. Therefore, in very specific cases the effect of He- 
enhancement may have to be considered in globular clusters. 

On the other hand for T e ff approaching the value of 3500 K 
and for most filters in the blue part of the spectrum (from U 
to V), \ABCs A \ become significantly larger and can amount 
to as much as 0.06 mag at [M/H] < -0.5, and 0.15 mag at 
[M/H] = +0.5. These low r e ff values are those typical of early- 
M giants, including for instance the tip of the RGB (TRGB) 
at old ages and moderately low metallicities ([Fe/H] ~ -0.7). 
Fortunately in the same range of low T e ff the red and infrared 
passbands present small ABCs corrections, in particular in the 
/ band ABC/ becomes smaller than 0.015 mag. We notice that 
distance determinations of resolved galaxies via the TRGB I- 
band magnitude should not be affected by possible galaxy-to- 
galaxy changes in the mean He content, since they usually refer 
to stars hotter than ~ 4000 K, for which the possible ABC/ cor- 
rections are even smaller than at 3500 K. 

However, we remark that it is not at all clear whether the 
significant values of | ABCs J at T s g ~ 3500 K are a serious 
problem owing to the well-known uncertainties of the ATLAS9 
models for < 4000 K. For instance, at T e $ ~ 3500 K starts 
the formation of strong molecular bands in the stellar spectra, 
which are not accurately reproduced - at least not at the level 
of a few percent - by present-day ATLAS9 models (see for 
instance Fluks et al. 1994). Among the reasons there is the lack 
in the line opacity computations of both triatomic molecules 
(with exception for H2O which is considered) and of numerous 
diatomic molecular transitions. 

Therefore, the significant changes in ABCs A that we find at 
low T e ff may be just one additional - and secondary - problem 
in a field that is already complicated in itself, and for which 
synthetic photometry has always been recognized not to pro- 
vide accurate answers. 

In conclusion, we find that the effects of changes in He 
abundances among stellar populations are quite modest when 
we look at the stellar atmospheres and their predicted bolo- 
metric corrections. Therefore, the use of tables of BCs com- 
puted for a single Y(Z) relation, is an acceptable approxima- 
tion in most cases. We provide tables for ABCs , in a series of 
[M/H], AY, and r e jf values, that may help the reader to evalu- 
ate whether this is an issue in the interpretation of their obser- 
vations. The effects of changing Y by as much as 0.1, instead, 
may have a quite high impact on the stellar evolutionary tracks, 
and have to be considered whenever it is suspected, as in the 
case of u> Cen. 
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Fig. 2. The quantity SF^ (see eq.[TJ as a function of A, as derived for the [M/H] = -0.5 spectra with AY - and AY = 0.1. The 
top and bottom panels show sequences of decreasing for both dwarfs and giants, respectively. Similar figures are available 
also for [M/H] = —1.5 and [M/H] = +0.5, and can be provided upon request. 
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Table 1. ABC Sa values (in mag) for the UBVRIJHK system of Bessell (1990) and Bessell & Brett (1988). 
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-0.008 


-0.010 


5000 


1.50 


0.001 


0.008 


0.008 


0.004 


0.004 


0.003 


0.000 


-0.009 


-0.012 


6000 


1.50 


-0.005 


0.002 


0.002 


0.001 


0.002 


0.001 


0.000 


-0.004 


-0.005 


8000 


1.50 


0.004 


-0.001 


-0.001 


-0.001 


0.000 


0.000 


0.000 


0.000 


0.000 










[M/H] 


= +0.5, AF = 0.2 










3500 


4.50 


0.060 


-0.107 


-0.103 


-0.147 


-0.102 


-0.058 


0.017 


0.042 


0.034 


4000 


4.50 


0.051 


-0.035 


-0.033 


-0.060 


-0.052 


-0.012 


0.009 


0.019 


0.019 


5000 


4.50 


0.013 


0.004 


0.003 


0.002 


-0.001 


-0.001 


0.000 


-0.002 


-0.002 


6000 


4.50 


0.001 


0.003 


0.003 


0.001 


0.000 


0.000 


0.000 


-0.002 


-0.003 


8000 


4.50 


-0.004 


0.006 


0.007 


0.004 


0.003 


0.001 


0.002 


0.000 


0.001 


12000 


4.50 


0.009 


0.005 


0.005 


0.000 


0.001 


0.001 


0.001 


0.002 


0.000 


20000 


4.50 


-0.002 


0.003 


0.003 


0.001 


0.000 


0.000 


-0.001 


-0.001 


-0.002 


40000 


4.50 


0.006 


0.004 


0.004 


0.007 


0.006 


0.004 


-0.002 


-0.003 


-0.003 


3500 


1.50 


-0.033 


-0.081 


-0.081 


-0.058 


-0.014 


-0.013 


0.006 


-0.002 


0.001 


4000 


1.50 


-0.012 


-0.006 


-0.006 


0.000 


0.002 


0.004 


0.006 


-0.007 


-0.007 


5000 


1.50 


-0.011 


0.005 


0.005 


0.004 


0.004 


0.003 


0.001 


-0.007 


-0.010 


6000 


1.50 


-0.011 


0.002 


0.001 


0.002 


0.002 


0.001 


0.000 


-0.003 


-0.004 


8000 


1.50 


0.010 


-0.001 


-0.001 


-0.001 


0.000 


0.000 


-0.001 


-0.001 


-0.001 



